PhytoREF: a reference database of the plastidial 16S rRNA gene of photosynthetic eukaryotes with curated taxonomy

Photosynthetic eukaryotes have a critical role as the main producers in most ecosystems of the biosphere. The ongoing environmental metabarcoding revolution opens the perspective for holistic ecosystems biological studies of these organisms, in particular the unicellular microalgae that often lack distinctive morphological characters and have complex life cycles. To interpret environmental sequences, metabarcoding necessarily relies on taxonomically curated databases containing reference sequences of the targeted gene (or barcode) from identified organisms. To date, no such reference framework exists for photosynthetic eukaryotes. In this study, we built the PhytoREF database that contains 6490 plastidial 16S rDNA reference sequences that originate from a large diversity of eukaryotes representing all known major photosynthetic lineages. We compiled 3333 amplicon sequences available from public databases and 879 sequences extracted from plastidial genomes, and generated 411 novel sequences from cultured marine microalgal strains belonging to different eukaryotic lineages. A total of 1867 environmental Sanger 16S rDNA sequences were also included in the database. Stringent quality filtering and a phylogeny-based taxonomic classification were applied for each 16S rDNA sequence. The database mainly focuses on marine microalgae, but sequences from land plants (representing half of the PhytoREF sequences) and freshwater taxa were also included to broaden the applicability of PhytoREF to different aquatic and terrestrial habitats. PhytoREF, accessible via a web interface (http://phytoref.fr), is a new resource in molecular ecology to foster the discovery, assessment and monitoring of the diversity of photosynthetic eukaryotes using high-throughput sequencing.

PhytoREF: a reference database of the plastidial 16S rRNA gene of photosynthetic eukaryotes with curated taxonomy

Photosynthetic eukaryotes have a critical role as the main producers in most ecosystems of the biosphere. The ongoing environmental metabarcoding revolution opens the perspective for holistic ecosystems biological studies of these organisms, in particular the unicellular microalgae that often lack distinctive morphological characters and have complex life cycles. To interpret environmental sequences, metabarcoding necessarily relies on taxonomically curated databases containing reference sequences of the targeted gene (or barcode) from identified organisms. To date, no such reference framework exists for photosynthetic eukaryotes. In this study, we built the PhytoREF database that contains 6490 plastidial 16S rDNA reference sequences that originate from a large diversity of eukaryotes representing all known major photosynthetic lineages. We compiled 3333 amplicon sequences available from public databases and 879 sequences extracted from plastidial genomes, and generated 411 novel sequences from cultured marine microalgal strains belonging to different eukaryotic lineages. A total of 1867 environmental Sanger 16S rDNA sequences were also included in the database. Stringent quality filtering and a phylogeny-based taxonomic classification were applied for each 16S rDNA sequence. The database mainly focuses on marine microalgae, but sequences from land plants (representing half of the PhytoREF sequences) and freshwater taxa were also included to broaden the applicability of PhytoREF to different aquatic and terrestrial habitats. PhytoREF, accessible via a web interface (http://phytoref.fr), is a new resource in molecular ecology to foster the discovery, assessment and monitoring of the diversity of photosynthetic eukaryotes using high-throughput sequencing.

PCR from the CPR offers a historical perspective on marine population ecology.

The Continuous Plankton Recorder (CPR) survey has collected plankton samples from regular tracks across the world's oceans for almost 70 y. Over 299,000 spatially extensive CPR samples are archived and stored in buffered formalin. This CPR archive offers huge potential to study changes in marine communities using molecular data from a period when marine pollution, exploitation and global anthropogenic impact were much less pronounced. However, to harness the amount of data available within the CPR archive fully, it is necessary to improve techniques of larval identification, to genus and species preferably, and to obtain genetic information for historical studies of population ecology. To increase the potential of the CPR database this paper describes the first extraction, amplification by the polymerase chain reaction and utilization of a DNA sequence (mitochondrial 16S rDNA) from a CPR sample, a formalin fixed larval sandeel.

Morphological and Phylogenetic Characterization of New Gephyrocapsa Isolates Suggests Introgressive Hybridization in the Emiliania/Gephyrocapsa Complex (Haptophyta)

The coccolithophore genus Gephyrocapsa contains a cosmopolitan assemblage of pelagic species, including the bloom-forming Gephyrocapsa oceanica, and is closely related to the emblematic coccolithophore Emiliania huxleyi within the Noëlaerhabdaceae. These two species have been extensively studied and are well represented in culture collections, whereas cultures of other species of this family are lacking. We report on three new strains of Gephyrocapsa isolated into culture from samples from the Chilean coastal upwelling zone using a novel flow cytometric single-cell sorting technique. The strains were characterized by morphological analysis using scanning electron microscopy and phylogenetic analysis of 6 genes (nuclear 18S and 28S rDNA, plastidial 16S and tufA, and mitochondrial cox1 and cox3 genes). Morphometric features of the coccoliths indicate that these isolates are distinct from G. oceanica and best correspond to G. muellerae. Surprisingly, both plastidial and mitochondrial gene phylogenies placed these strains within the E. huxleyi clade and well separated from G. oceanica isolates, making Emiliania appear polyphyletic. The only nuclear sequence difference, 1 bp in the 28S rDNA region, also grouped E. huxleyi with the new Gephyrocapsa isolates and apart from G. oceanica. Specifically, the G. muellerae morphotype strains clustered with the mitochondrial β clade of E. huxleyi, which, like G. muellerae, has been associated with cold (temperate and sub-polar) waters. Among putative evolutionary scenarios that could explain these results we discuss the possibility that E. huxleyi is not a valid taxonomic unit, or, alternatively the possibility of past hybridization and introgression between each E. huxleyi clade and older Gephyrocapsa clades. In either case, the results support the transfer of Emiliania to Gephyrocapsa. These results have important implications for relating morphological species concepts to ecological and evolutionary units of diversity.

Morphological and Phylogenetic Characterization of New Gephyrocapsa Isolates Suggests Introgressive Hybridization in the Emiliania/Gephyrocapsa Complex (Haptophyta)

The coccolithophore genus Gephyrocapsa contains a cosmopolitan assemblage of pelagic species, including the bloom-forming Gephyrocapsa oceanica, and is closely related to the emblematic coccolithophore Emiliania huxleyi within the Noëlaerhabdaceae. These two species have been extensively studied and are well represented in culture collections, whereas cultures of other species of this family are lacking. We report on three new strains of Gephyrocapsa isolated into culture from samples from the Chilean coastal upwelling zone using a novel flow cytometric single-cell sorting technique. The strains were characterized by morphological analysis using scanning electron microscopy and phylogenetic analysis of 6 genes (nuclear 18S and 28S rDNA, plastidial 16S and tufA, and mitochondrial cox1 and cox3 genes). Morphometric features of the coccoliths indicate that these isolates are distinct from G. oceanica and best correspond to G. muellerae. Surprisingly, both plastidial and mitochondrial gene phylogenies placed these strains within the E. huxleyi clade and well separated from G. oceanica isolates, making Emiliania appear polyphyletic. The only nuclear sequence difference, 1 bp in the 28S rDNA region, also grouped E. huxleyi with the new Gephyrocapsa isolates and apart from G. oceanica. Specifically, the G. muellerae morphotype strains clustered with the mitochondrial β clade of E. huxleyi, which, like G. muellerae, has been associated with cold (temperate and sub-polar) waters. Among putative evolutionary scenarios that could explain these results we discuss the possibility that E. huxleyi is not a valid taxonomic unit, or, alternatively the possibility of past hybridization and introgression between each E. huxleyi clade and older Gephyrocapsa clades. In either case, the results support the transfer of Emiliania to Gephyrocapsa. These results have important implications for relating morphological species concepts to ecological and evolutionary units of diversity.

Unexpected prevalence of parasite 18S rDNA sequences in winter among Antarctic marine protists

Parasites are not typically considered to be important components of polar marine ecosystems. It was therefore surprising when 18S rDNA surveys of protists in the West Antarctic Peninsula in winter revealed high abundances of parasite sequences. Parasite sequences made up, on average, over half (52%) of sequence reads in samples from deep water in winter. Winter surface water and sediment samples contained relatively fewer, but still strikingly high, parasite sequence reads (13 and 9%, respectively), while surface water samples in summer contained fewer parasite sequences (1.8%). A total of 1028 distinct parasite Operational Taxonomic Units were observed in winter, with the largest abundances and diversities within Syndiniales groups I and II, including Amoebophrya. Less abundant parasite sequence groups included Apicomplexa, Blastodinium, Chytriodinium, Cryptocaryon, Paradinium, Perkinsidae, Pirsonia and Ichthyophonae. Parasite sequence distributions suggested interactions with known hosts, such as diatom parasites which were mainly in the sediments, where resting spores of Chaetoceros spp. diatoms were abundant. Syndiniales sequences were correlated with radiolarian sequences, suggesting parasite–host interactions. The abundant proportions of parasite sequences indicate a potentially important role for parasites in the Antarctic marine ecosystem, with implications for plankton population dynamics, the role of the microbial loop, carbon flows and ecosystem responses to ongoing anthropogenic climate change.

Unexpected prevalence of parasite 18S rDNA sequences in winter among Antarctic marine protists

Parasites are not typically considered to be important components of polar marine ecosystems. It was therefore surprising when 18S rDNA surveys of protists in the West Antarctic Peninsula in winter revealed high abundances of parasite sequences. Parasite sequences made up, on average, over half (52%) of sequence reads in samples from deep water in winter. Winter surface water and sediment samples contained relatively fewer, but still strikingly high, parasite sequence reads (13 and 9%, respectively), while surface water samples in summer contained fewer parasite sequences (1.8%). A total of 1028 distinct parasite Operational Taxonomic Units were observed in winter, with the largest abundances and diversities within Syndiniales groups I and II, including Amoebophrya. Less abundant parasite sequence groups included Apicomplexa, Blastodinium, Chytriodinium, Cryptocaryon, Paradinium, Perkinsidae, Pirsonia and Ichthyophonae. Parasite sequence distributions suggested interactions with known hosts, such as diatom parasites which were mainly in the sediments, where resting spores of Chaetoceros spp. diatoms were abundant. Syndiniales sequences were correlated with radiolarian sequences, suggesting parasite–host interactions. The abundant proportions of parasite sequences indicate a potentially important role for parasites in the Antarctic marine ecosystem, with implications for plankton population dynamics, the role of the microbial loop, carbon flows and ecosystem responses to ongoing anthropogenic climate change.

Defining seasonal marine microbial community dynamics

Here we describe, the longest microbial time-series analyzed to date using high-resolution 16S rRNA tag pyrosequencing of samples taken monthly over 6 years at a temperate marine coastal site off Plymouth, UK. Data treatment effected the estimation of community richness over a 6-year period, whereby 8794 operational taxonomic units (OTUs) were identified using single-linkage preclustering and 21 130 OTUs were identified by denoising the data. The Alphaproteobacteria were the most abundant Class, and the most frequently recorded OTUs were members of the Rickettsiales (SAR 11) and Rhodobacteriales. This near-surface ocean bacterial community showed strong repeatable seasonal patterns, which were defined by winter peaks in diversity across all years. Environmental variables explained far more variation in seasonally predictable bacteria than did data on protists or metazoan biomass. Change in day length alone explains >65% of the variance in community diversity. The results suggested that seasonal changes in environmental variables are more important than trophic interactions. Interestingly, microbial association network analysis showed that correlations in abundance were stronger within bacterial taxa rather than between bacteria and eukaryotes, or between bacteria and environmental variables.

Ocean acidification shows negligible impacts on high-latitude bacterial community structure in coastal pelagic mesocosms

The impact of ocean acidification and carbonation on microbial community structure was assessed during a large-scale in situ costal pelagic mesocosm study, included as part of the EPOCA 2010 Arctic campaign. The mesocosm experiment included ambient conditions (fjord) and nine mesocosms with pCO(2) levels ranging from similar to 145 to similar to 1420 mu atm. Samples for the present study were collected at ten time points (t-1, t1, t5, t7, t12, t14, t18, t22, t26 to t28) in seven treatments (ambient fjord (similar to 145), 2x similar to 185, similar to 270, similar to 685, similar to 820, similar to 1050 mu atm) and were analysed for "small" and "large" size fraction microbial community composition using 16S rRNA (ribosomal ribonucleic acid) amplicon sequencing. This high-throughput sequencing analysis produced similar to 20 000 000 16S rRNA V4 reads, which comprised 7000OTUs. The main variables structuring these communities were sample origins (fjord or mesocosms) and the community size fraction (small or large size fraction). The community was significantly different between the unenclosed fjord water and enclosed mesocosms (both control and elevated CO2 treatments) after nutrients were added to the mesocosms, suggesting that the addition of nutrients is the primary driver of the change in mesocosm community structure. The relative importance of each structuring variable depended greatly on the time at which the community was sampled in relation to the phytoplankton bloom. The sampling strategy of separating the small and large size fraction was the second most important factor for community structure. When the small and large size fraction bacteria were analysed separately at different time points, the only taxon pCO(2) was found to significantly affect were the Gammaproteobacteria after nutrient addition. Finally, pCO(2) treatment was found to be significantly correlated (non-linear) with 15 rare taxa, most of which increased in abundance with higher CO2.

Minor impact of ocean acidification to the composition of the active microbial community in an Arctic sediment

Effects of ocean acidification on the composition of the active bacterial and archaeal community within Arctic surface sediment was analysed in detail using 16S rRNA 454 pyrosequencing. Intact sediment cores were collected and exposed to one of five different pCO(2) concentrations [380 (present day), 540, 750, 1120 and 3000 atm] and RNA extracted after a period of 14 days exposure. Measurements of diversity and multivariate similarity indicated very little difference between pCO(2) treatments. Only when the highest and lowest pCO(2) treatments were compared were significant differences evident, namely increases in the abundance of operational taxonomic units most closely related to the Halobacteria and differences to the presence/absence structure of the Planctomycetes. The relative abundance of members of the classes Planctomycetacia and Nitrospira increased with increasing pCO(2) concentration, indicating that these groups may be able to take advantage of changing pH or pCO(2) conditions. The modest response of the active microbial communities associated with these sediments may be due to the low and fluctuating pore-water pH already experienced by sediment microbes, a result of the pH buffering capacity of marine sediments, or due to currently unknown factors. Further research is required to fully understand the impact of elevated CO2 on sediment physicochemical parameters, biogeochemistry and microbial community dynamics.

Morphological and genetic variation in the north Atlantic copepod, Centropages typicus

This study describes phenotypic and genotypic variations in the planktonic copepod, Centropages typicus (Copepoda: Calanoida) that indicate differentiation between geographical samples. We found consistent differences in the morphology of the chela of the sexually modified fifth pereiopod (P5) of male C. typicus between samples from the Mediterranean, western North Atlantic and eastern North Atlantic. A 560 base pairs (bp) region of the C. typicus mitochondrial cytochrome c oxidase subunit I (COI) and a 462 bp fragment of the nuclear rDNA internal transcribed spacer (ITS) tandem array were analysed to determine whether these morphological variations reflect population genetic differentiation. Mitochondrial haplotype diversity was found to be high with 100 unique COI haplotypes among 116 individuals. Analysis of mtCOI variation suggested differentiation between the Mediterranean and Atlantic populations but no separation was detected within the Atlantic. Intragenomic variation in the ITS array suggested genetic differentiation between samples from the western North Atlantic and those from the eastern North Atlantic and Mediterranean. Breeding experiments would be required to elucidate the extent of genetic isolation between C. typicus from the different population centres.

Diverse profiles of N‐acyl‐homoserine lactone molecules found in cnidarians

Many marine habitats, such as the surface and tissues of marine invertebrates, including corals, harbour diverse populations of microorganisms, which are thought to play a role in the health of their hosts and influence mutualistic and competitive interactions. Investigating the presence and stability of quorum sensing (QS) in these ecosystems may shed light on the roles and control of these bacterial communities. Samples of 13 cnidarian species were screened for the presence and diversity of N-acyl-homoserine lactones (AHLs; a prevalent type of QS molecule) using thin-layer chromatography and an Agrobacterium tumefaciens NTL4 biosensor. Ten of 13 were found to harbour species-specific, conserved AHL profiles. AHLs were confirmed in Anemonia viridis using liquid chromatography tandem mass spectrometry. To assess temporal role and stability, AHLs were investigated in A. viridis from intertidal pools over 16 h. Patterns of AHLs showed conserved profiles except for two mid-chain length AHLs, which increased significantly over the day, peaking at 20:00, but had no correlation with pool chemistry. Denaturing gel electrophoresis of RT-PCR-amplified bacterial 16S rRNA showed the presence of an active bacterial community that changed in composition alongside AHL profiles and contained a number of bands that affiliate with known AHL-producing bacteria. Investigations into the quorum sensing-controlled, species-specific roles of these bacterial communities and how these regulatory circuits are influenced by the coral host and members of the bacterial community are imperative to expand our knowledge of these interactions with respect to the maintenance of coral health.

Response of an Arctic Sediment Nitrogen Cycling Community to Increased CO2

Ocean acidification influences sediment/water nitrogen fluxes, possibly by impacting on the microbial process of ammonia oxidation. To investigate this further, undisturbed sediment cores collected from Ny Alesund harbour (Svalbard) were incubated with seawater adjusted to CO2 concentrations of 380, 540, 760, 1,120 and 3,000 μatm. DNA and RNA were extracted from the sediment surface after 14 days' exposure and the abundance of bacterial and archaeal ammonia oxidising (amoA) genes and transcripts quantified using quantitative polymerase chain reaction. While there was no change to the abundance of bacterial amoA genes, an increase to 760 μatm pCO2 reduced the abundance of bacterial amoA transcripts by 65 %, and this was accompanied by a shift in the composition of the active community. In contrast, archaeal amoA gene and transcript abundance both doubled at 3,000 μatm, with an increase in species richness also apparent. This suggests that ammonia oxidising bacteria and archaea in marine sediments have different pH optima, and the impact of elevated CO2 on N cycling may be dependent on the relative abundances of these two major microbial groups. Further evidence of a shift in the balance of key N cycling groups was also evident: the abundance of nirS-type denitrifier transcripts decreased alongside bacterial amoA transcripts, indicating that NO3 − produced by bacterial nitrification fuelled denitrification. An increase in the abundance of Planctomycete-specific 16S rRNA, the vastmajority of which grouped with known anammox bacteria, was also apparent at 3,000 μatm pCO2. This could indicate a possible shift from coupled nitrification–denitrification to anammox activity at elevated CO2.

Live discrimination of Calanus glacialis and C. finmarchicus females: can we trust phenological differences?

Two key players in the Arctic and subarctic marine ecosystem are the calanoid copepods, Calanus finmarchicus and C. glacialis. Although morphologically very similar, these sibling species have different life cycles and roles in the Arctic pelagic marine ecosystem. Considering that the distribution of C. glacialis corresponds to Arctic water masses and C. finmarchicus to Atlantic water masses, the species are frequently used as climate indicators. Consequently, correct identification of the two species is essential if we want to understand climate-impacted changes on Calanus-dominated marine ecosystems such as the Arctic. Here, we present a novel morphological character (redness) to distinguish live females of C. glacialis and C. finmarchicus and compare it to morphological (prosome length) and genetic identification. The characters are tested on 300 live females of C. glacialis and C. finmarchicus from Disko Bay, western Greenland. Our analysis confirms that length cannot be used as a stand-alone criterion for separation. The results based on the new morphological character were verified genetically using a single mitochondrial marker (16S) and nuclear loci (six microsatellites and 12 InDels). The pigmentation criterion was also used on individuals (n = 89) from Young Sound fjord, northeast Greenland to determine whether the technique was viable in different geographical locations. Genetic markers based on mitochondrial and nuclear loci were corroborative in their identification of individuals and revealed no hybrids. Molecular identification confirmed that live females of the two species from Greenlandic waters, both East and West, can easily be separated by the red pigmentation of the antenna and somites of C. glacialis in contrast to the pale opaque antenna and somites of C. finmarchicus, confirming that the pigmentation criterion is valid for separation of the two species